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In this work, we have investigated whether Y(4260) and Z J (4250) could be DiD or DoD* molecules 
in the framework of meson exchange model. The off-diagonal interaction induced by n exchange 
plays a dominant role. The o exchange has been taken into account, which leads to diagonal 
interaction. The contribution of a exchange is not favorable to the formation of molecular state 
^ ■ with I'^(J'''^) = 0"(1 ), however, it is beneficial to the binding of molecule with = 1" (1"). 

I Light vector meson exchange leads to diagonal interaction as well. For Z^(4250), the contribution 

from p and u exchange almost cancels each other. For the currently allowed values of the effective 
coupling constants and a reasonable cutoff A in the range 1-2 GeV, We find that Y(4260) could be 
^ ' accommodated as a DiD and DoD* molecule, whereas the interpretation of Z J (4250) as a DiD or 

O ■ DqD* molecule is disfavored. The bottom analog of Y(4260) and Z J (4250) may exist, and the most 

' promising channels to discovery them are 7r^7r~T and iv^Xbi respectively. 
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PACS numbers: 12.39.Pn, 12.40.Yx, 13.75.Lb,12.39.Hg 



^ . I. INTRODUCTION 

^ i' In the past years, a number of charmonium-like X, Y, Z states have been observed, which stimulate a lot of discussion 
D \ about the structures and properties of these resonances. In particular, the Z+(4430) observed in the ■K'^ij}' invariant 
/ spectrum[l| carries one unit electric charge. Consequently, it can not be simple charmonium. Recently, two new 

resonance-like structures Z^(4051) and Zo'(4250) in the tt~^Xc1 mass distribution in exclusive B — > K~n~^Xci have 
' been reported by the Belle collaboration [J- Their masses and widths arc determined to be Mi = (4051 ± Utli)McV, 
^ ' Ti = (82j:^^t^^)MeV, Ms (4248l29i35°)McV and T2 = (177;^gt^i'^)McV respectively, with the product branching 
fractions K- Z+^) x S(Z+2 ^ tt+Xci) = (S.OtJjt?;^) x 10"^ and {^..^tUt'^o^ ) x 10"^ respectively. Both 

Z^(4051) and Z^(4250) carry one unit electric charge like Z+(4430), hence they must be states beyond quark model, 
if these states are confirmed in future. Since tt"'' is an isovector with negative G-parity, and Xci is a isospin singlet 
0^ ■ with positive G-parity, the quantum numbers of Z^(4051) and Zj(4250) are I*^ = 1^. It is remarkable that some 
, states are in the vicinity of the S-wave threshold of two charmed mesons, e.g., X(3872) and Z+(4430) are very 
close to the thresholds of D*D and DiD* respectively, therefore it is tempting to interpret these states as molecular 
statesjl, II] . Particularly, Y(4260) and Z^(4250) are close to the DiD and DqD* thresholds, which inspires the 
theoretical interpretations of Y(4260) as a DqD* molecule and Z+(4250) as a DiD molecule[6|. 

Y(4260) was reported by the Babar collaboration in the -k^tx^ J ji]} invariant spectrum of the reaction e^e~ 
' Iisrt^^t^' J li'^ii which has been confirmed by both the CLEO and Belle collaboration[^Q. A fit to the peak with 
■ a single Breit- Winger resonance shape yields a mass M — (4259 ± 10)McV and the full width F = (88 ± 24)McV. 
Evidently the state is a vector with cc flavor, and its quantum numbers are determined to be l'-'(J^'^) = 0^(1^^). 
Although it is above the threshold for decaying into DD, DD*(D*D) or D*D* meson pairs, there is no evidence for 
Y(4260) in these channels [H [ul, [12. Therefore Y(4260) appears not to be a canonical charmonium. 

The observation of the Y(4260) has sparked many theoretical speculations. It has variously been identified as a 
conventional V'(4'S') based on a relativistic quark model[l3l|. a tetraquark cess state[l3| which decays predominantly 
into DsDg, or a charmonium hybrid [isj. The data on e+e" — *■ DgDs show a peaking above threshold around 4 GeV 
but no evidence of afhnity for a structure at 4.26GeV[l^. If these data are confirmed, then the interpretation of 
Y(4260) as a cess tetraquark would be ruled out. Moreover, dynamical calculation of tetraquark states indicated that 
Y(4260) can not be interpreted as P-wave 1 state of charm-strange diquark-antidiquark, because the corresponding 
mass is found to be 200 MeV heavier p7j. Although the charmonium hybrid is a very attractive interpretation, the 
Lattice QCD simulations predict that lightest charmonium hybrid is about 4.4GeV[i3|, which is very close to the 
new charmonium-like state Y(4360)[l5. As has been proposed in Ref.jl^l, a possible resolution to this issue is that 
Y(4360) is the candidate of charmonium hybrid, while Y(4260) is a DiD hadronic molecule. 

In Ref . [2]| . Swanson emphasized that we should examine the DiD molecular interpretation before finally concluding 
that Y(4260) is a charmonium hybrid. Furthermore, he pointed out that tt exchange does not lead to a diagonal 
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interaction in the DiD channel, and certain novel mechanism such as ofF-diagonal interaction may be required. In 
Refs.(22l. [23j . Close showed that parity conservation requires the tt vertex to link D D* and Di Dg, then the tt 
exchange gives an off-diagonal potential linking DiD DqD* or DDi ^ D*Do. This tt exchange attraction possibly 
results in a 1~~ hadronic molecule near the DiD threshold. In this work, we shall investigate whether Y(4260) and 
Z J (4250) could be hadronic molecule due to the off-diagonal tt exchange effect in the framework of heavy quark 
effective theory. The contribution of a exchange has been considered, which results in diagonal interaction. The light 
vector mesons p and uj exchange is discussed as well. 

The paper is organized as follows. In section II, we present the formalism to include both heavy meson and 
anti- meson fields in the heavy meson chiral perturbation theory(HMxPT), and the complete Lagrangian is written 
out explicitly. Section III illustrates the systematic procedure for converting a general T-matrix into an equivalent 
potential operator. Later wc follow this to derive the effective potential. In section IV, we present both the diagonal 
and non-diagonal potential related with Y(4260) and Z J (4250). In section V, we investigate the possible bound states 
of the DiD and DpD* system by solving the coupled-channel Schrodinger equations, and the structures of Y(4260) 
and Z^(4250) are discussed. Moreover, the bottom analog of Y(4260) and Z J (4250) is studied. We present our 
conclusions and some relevant discussions in Section VI. Finally, the potential from p and cu exchange is shown in the 
Appendix. 



II. FORMALISM FOR THE SYSTEM CONTAINING BOTH MESON AND ANTI-MESON FIELDS IN 

HMxPT 



The strong interaction between pseudo-Goldstone bosons and the mesons containing a heavy quark is described 
by the so-called heavy meson chiral perturbation theory (HMxPT) [H, [2^, . The heavy meson chiral perturbation 
theory is constructed starting from the spin-flavor symmetry occurring in QCD in the infinite heavy quark mass limit, 
and from the chiral symmetry valid in the massless limit for the light quarks. In HM^PT, the heavy-light meson field 
appears in a covariant form, which is represented by a 4 x 4 Dirac-type matrix. The negative and positive parity 
doublets containing a heavy quark Q and a light anti-quark of flavor a, can be respectively described by the superfields 
Ha, Sa and as follows 



(1) 



The above various operators annihilate mesons of four-velocity v which is conserved in strong interaction pro- 
cesses. The heavy field operators contain a factor v^Mp and have dimension 3/2. Under a heavy quark spin SU{2) 
transformation S and a generic light flavor transformation U (i.e., U £ SU{3)) 



HiQ) 



(2) 



The conjugate field, which creates heavy-light mesons containing a heavy quark Q and a light anti-quark of flavor a, 
is defined as 



which transforms under 5* and U as 

TjiQ') TT TjiQ') ci -c(Q) TT -niQ) at 7r(Q)f^ tt 7f;(Q)Mr.t 
-"a ^ Uabilb "J : a —^^ab^b i -'a —* ^ ab^ b 



(3) 



(4) 



The octet of light pseudoscalar mesons can be introduced using the non-linear representation E = and ^ = 
exp(iA^//7r) with /^r = 132 MeV. The matrix M contains tt, K, rj fields, which is a 3 x 3 hermitian and traceless 
matrix 



M 



K- 



.11 + Jl. K" 
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Under the chiral symmetry, the field ^ transforms as 

e ^ .9LCC/t = (7Cg]j (6) 

where and p^j are left-handed and right-handed global SU{3) transformation respectively. 

The effective QCD Lagrangian is constructed b y im posing invariance under both heavy quark spin-flavor transfor- 
mation and chiral transformation, it is [2^. [25l. [26l. l27l. 

+zMT(«'^A.,.75^i''^) + *^(rf ^''(i^,4)fea75i?f'=^) + z^(rf )''(^^,),„75i?i'^^) + /^.c.] (7) 

where (...) means trace over the 4x4 matrices, the covariant derivative = 9^ + V^, the vector current and the 
axial current Af^ are defined by 

■A,^lied,^-Cd,i^) (8) 

In order to describe mesons containing heavy anti-quark Q, we have to introduce six new fields Pal!^\ Pa^\ Pifj} i 
-^Oa I ^20^1/ ^^"i ^la/i which dcstroy mesons containing a heavy anti-quark Q. The phase of the field Pa^ relative 
to Pal?^ , Pa^^ to Pi*^^ etc can be fixed by the following charge conjugation convention 

p*(Q) _ rp*{Q)r-^ p*(Q) _ _r p*(Q)r~^ p(Q) _ n p'^^^r-'^ 

CiC-^^e, CV,,C-^ = -Vl, CA,,C-'^Al (9) 



The mesons containing a heavy anti-quark Q and a light quark of flavor a can be included into the theory by applying 

srfieldi 

1- 



the charge conjugation operation to the above heavy-light meson superfields Hli^\ si'^^ and 
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Ti«- ^ C(Cri^)C-i)^C-i [Pi^)'^^. - y^PL^j75(5'^^ - i(7'^ - «^)7^)]^^ (10) 

The matrix C is the charge conjugation matrix for Dirac spinors with C = 47^7'^, and the transpose is on the spinor 
matrix indices. Under the heavy quark spin transformation S and light quark SU{3) flavor symmetry U, 

Hp-.Ua,HpS\ Sp^UabSpS\ T(Q^-*Uai,Tg^S^ (11) 
Similarly the hermitian conjugate flelds are defined by 

H'^'^^oHpho, Sp=ioSpho, ri?=7oTi?t^„ (12) 
Under the symmetry transformation 5* and U 

Hp^SHpul Sp^SSpul rS'^STifuL (13) 

For the system including both heavy meson and heavy anti- meson field in the HMxPT, the total effective Lagrangian 
should be invariant under the charge conjugation transformation. The interaction between the pseudo-Goldstone 
bosons and the meson containing one heavy anti-quark can be obtained from Eq.([7]) by applying the charge conjugation 
operator 



+zh{SpA,a,l,Tpn + ^^{H^?\4 + ^^{H?\a, ^')a,j,Tp') + h.c. (14) 
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where D'^ = d^~V^. After expanding the effective Lagrangian in Eq.® and Eq. (fT4)) to the leading order of pseudo- 
Goldstone field, we can obtain the following effective interactions, which is needed in our work 

^DoDiP .9DoDipD^b(a^Al)baDj„ + .gp^D^p V)la{d^,M)abGlb + h.C. 
>CdDoP = i5DDoP(Doh 9m ^i)dt,Mba + ».9dDop(^06 9m ^\)d^Mab + /i-c 
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Md-Mdi 



d^T^'^,,{d,drM)bad^'D*l 



Md.Mdi 

In the chiral and heavy quark limit, the above coupling constants are 

2.9 



h.C. 



5dd'P - -.9dd*p 

ffDoDiP = 5doDiP = 
ffDDoP = 5dDoP = 



2^6 h 



VMdMe 



3 U 



^MdoMdi 



a/6 hi + /i2 



ffD-DiP = ffD*DiP = 3- ^ J a/Md'Mdi 



(15) 



(16) 



We would like to stress that the DD*P coupling constant is the negative of the DD P coupling constant, because 
of the phase convention for charge conjugation chosen in Eq.®. The effective Lagrangian between a and heavy 
nieson(anti-nieson) are [13] 



h' 



Jtt Jtt 



The coupling constants are estimated as follows j3C 

.97r 



2\/6' 



(17) 



(18) 



where g^^ = 3.73 and gA = 0.6. As in Ref . [31| . we take \g'^\ ~ \ga-\ and = \ha-\ approximately when performing the 
numerical analysis. Expanding the Lagrangian £0-1 we get the interactions associated with a 



(19) 
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The relevant coupling constants are 



5DDcr — 


■9d D CT ~ 






yDiDi cr 


= -2(7"Mn 








5DDict = 


~ -^D Di cr = 
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~ 5D*D*cr 




ffDoDoCT 


= -^DoDoCT 


= 2.g;MDo 








ffD*Do(T 




(T f ^ 



(20) 



III. CONVERTING THE T-MATRIX INTO THE EFFECTIVE POTENTIAL 

The T-matrix for ^(pi)i?(p2) C(p']^)-D(p2) scattering process can be represented by an equivalent Born-order 
potential operator Vf,„(ri — r2, Vi, V2) between pointlike particles, the definition of this potential operator is [ssl. [36j 

^^(P'l +P2 - Pi -P2)T/,(Pl,P2,Pl,P2) 

1 



(27r)3 



cPrid^r2e-'^P'^ ''^+P'^ ''^^Vbn{ri - r2, Vi, V2)e*^Pi ""^+P= ''=) (21) 



where T/i(pi, p2, p'l, P2) is the T-matrix for the process A{px)B{p2) C(p'i)D(p2)- In general, Tji depends on 
all the involved momentum pi, p2, p'^ and P2. For convenience, we introduce 

Pi = ^(pi+p'l), P2 = ^(P2+P2), q = p'i Pl = P2 P2 (22) 

In the center of mass frame Pi = — P2- The amplitude Tfi can be expanded as a power series in Pi^ and P2i 

T/.(pi,P2,pi,P^) = T(")(q) +Tf^"'(q)Pn +Tf^)(q)P2, -f Tg^^)(q)PHP2, + ... (23) 

This procedure produces the full Breit-Fermi Hamiltonian when it is applied to the photon exchanged electron-electron 
scattering amplitude expanded to O(P^). The leading term T("'(q) is a function of q only, its Fourier transformation 
gives us a local potential V{r) that is a function of ri — r2 = r only. The relation between T'"' (q) and V{r) is 

^^'■^ = (2^)^/ ^24) 

For the higher terms of the T-matrix expansion. Pi,, and P2i are replaced by left- and right-gradients in the equivalent 
potential operator defined implicitly by Eg . (PT|) [SSj . Following this systematic procedure, we can convert a general 
T-matrix into an equivalent potential operator. In this work, we obtain the local potential by Fourier transforming 
the leading terms T^°)(q) of the scattering amplitude Tfi, which is common in potential model [37i. [ssj. 

Since the propagators are off-shell, we introduce form factor at each vertex when writing out the scattering ampli- 
tude, the usual form factor is expressed asfsol. [40j 

- ^ <-) 

where A is an adjustable constant within a reasonable range of 1-2 GeV, which models the off-shell effects at the 
vertices due to the internal structure of the meson, m and q are the mass and the four momentum of the exchanged 
meson respectively. 



IV. THE EFFECTIVE POTENTIALS RELATED WITH Y(4260) AND Z+(4250) 

Recently, the meson exchange model based on the HM^PT has been used to study possible heavy flavor molecule [Sll. 
l4l| . In this section, we will follow the general procedure shown above to derive the effective potential associated with 
Y(4260) and Z+(4250) in the framework of HMxPT. Fr om the effective interaction in Eg. p^ and Eq. p^ . we can 
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write down the corresponding scattering amplitude for each diagram, including the form factor at each vertex. Then 
we get the equivalent potential in momentum space following the general formalism presented in section III. Finally we 
make Fourier transformation to derive the potentials in coordinate space. Because of parity conservation, pseudoscalar 
TT and 77 exchange only contributes to the off-diagonal interaction, whereas a exchange and light vector mesons p, uj 
exchange result in diagonal interaction only. The corresponding scattering diagrams are shown in FiglT] 



D*(D ) Do(Do) Do(Do) D (D*) 



7r,?7 / \ 7r,77 



D(D) Di(Di) D(D) Di(Di) 

(a) {h) 



D(D) Di(Di) Di(Di) D(D) 



D(D) Di(Di) D(D) Di(Di) 

(c) (d) 



D*(D ) Do(Do) Do(Do) D (D*) 

a,p,uj/ \ a,p,uj 



D*(D*) Do(Do) D*(D*) Do(Do) 

(e) (/) 



FIG. 1: The scattering diagrams with pseudoscalars tt, 77 exchange, a exchange, and light vector mesons p, to exchange. 



Under the ansatz of Y(4260) as a DiD or DoD* hadronic molecule, we can write down its flavor wavefunction 

|Y(4260)) = i[|D;D°) + ID+D-) - |D°Di) - |D+Dr)] 

|Y'(4260)) = i[|Dj]D*") + |D+D*-) + |D*"Do") + ID^+D^-)] (26) 

We stress that the phase convention imder charge conjugation is consistent with Eq.Q. In the same way, the flavor 
wavefunction of Z^(4250) is 

|Z+(4250)) = i=[|D+D°) + |D+D°)] 

|Z^+(4250)) = ^[|Do+D*°) - |D*+D°)] (27) 
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In this case, its quantum number are I^(J^) = Following the procedure discussed above, we can calculate 

the effective potential in momentum space, it is a lengthy and tedious calculation. 
For Y(4260); the exchange potential in momentum space is 



q2 + A2 y q2 + ml 9/2 V q2 + / q2 + ^2 

-^^) ^rU^ihwT^^) ^T7i ^^^^ 

where we have include the monopole form factor in Eg. ipS]) to regularize the potential. The diagonal potential Fii(q) 
and V22(q) is induced by a exchange, and the non-diagonal potential Vi2(q) ( or V2i(q) ) arises from the pseudo- 
Goldstone bosons tt and 77 exchange. The effective potential from p, lo exchange is shown in the Appendix. The 
potential for Z J (4250) in momentum space is 



^ 76gfe/ A2-m2 N2 q2 ^6 g/i / A' - ^2 x 2 q2 



^l2(q) V2M) i8/2Vq2 + x2; q2+/,2 ' 54/2Vq2+X2; q2 + ^2 

'A2-m2x2 3,5^' 2/i;2.A2-to2 2 q2 



V ( } ^ ( '' Y ^"^'^ ( 

'^llW Vq2+A2y' q2+m2 9/2 V q2 + ATf 7 q2 + ^2 

/A2-m2x2 g^g' hi ( - ml\^ q2 , ^ 

^ ' V q2 + A2 / q2 + to2 3/2 V q2 + ATf / q2 + /il > 

The various parameters appearing in the above formulas are defined as follows. 
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-Mdo) 


f4 = 










^ll = 




-(Mdo 







(30) 

After performing Fourier transformation, we obtain the potential forms in configuration space, For Y(4260), the 
potential in coordinate space is 

Vi2{r) = V2i{r) = — — Z{k,Xi,^ii,m^,r) + — ■—Z{K,Xi,^i2,'m^,r) 



2h 
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Vii (r) = ga9aH{A, m,,r) + ^^(A, X2, Ms, m,, r) 



_?l2 

The potential in coordinate space for Z J (4250) is 

--Z(A,Xi,Mi,m.,r) + — - 

2^'2 
9/2 

_?l2 

3/.2 

Here the functions H{A, to, r) and Z{A, X, fi, to, r) are defined as 



V22(?') = gag'crH{A, m^r, r) + ^Z(A, X3, fi^, to^-, r) (31) 



Vi2ir) = V2i{r) = - — |^Z(A, Xi, /n, to,, r) + — ^Z(A, Xi, ^2, m^, r) 
Vii(r) = g„g'^H{A,ma,r) - ^^Z(A, X2, /is, to^, r) 

/l2 

^22(0 - .9^.9;i/(A,m„r)- ^Z(A,X3,/i4,m„,r) (32) 



i7(A, m, r) = ^^(e-"'- - e"^'^) - ^\ "'' g-^'^ 
^ ^ 47rr^ ^ SttA 

Z(A,X,/.,m,r) = ]--{X'e-''^- - f^'e-^ + ^^^^{X ^)e-^'^ (33) 
Air r Htt r 
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We take the typical values of the coupling constants gh = 0.85, gaQa = 0.58, gag'l ~ 0.58 and \ha\ = \h'^\ = 0.35, 
and A = 1.5GeV is chosen for an illustration, the variation of the effective potential with respect to r is shown in 
FiglSl It is obvious that the magnitude of the diagonal potentials from a exchange is smaller than that of the off- 
diagonal potential from tt and rj exchange, this is mainly because is small than m(j. Moreover, the magnitude 
of the off-diagonal potential related with Y(4260) are larger that associated with Z^(4250), the latter is about one 
third of the former. This is consistent with results from chiral quark model[3^, consequently the I'^(J^'^) = 0^(1^^) 
configuration is easier to bind than the = l-(l^) confi guration. 




1 2 3 4 1. 1.5 2.0 2.5 3.0 3.5 4.0 

r(fm) i-(fm) 



(c) (d) 

FIG. 2: (color onIine)The effective potential for the DiD and DqD* system from pseudoscalar tt, rj exchange and scalar a 
exchange. Solid line represents the non-diagonal potential Vi2(r) (or V2i(r)), short dashed and dash dotted lines respectively 
correspond to the diagonal potential Vii(j') and V22{r). (a) and (b) are related with Y(4260), and (b) shows the long range 
behavior of the potential, (c) and (d) are related with Z J (4250), and (d) is the long rang shape of the potential. 



V. THE STRUCTURES OF Y(4260) AND Z J (4250) AND THE BOTTOM ANALOG 

A. The bound states of the DiD and DqD* system with the structure of Y(4260) and Z J (4250) 

With the above effective potential, we shall explore whether there are bound states with I^(J^^) = 0^(1 ) or 
I°(JP) ^ 1^(1") in the DiD and DqD* system, by means of solving the two channels coupled Schrodinger equation. 
There are various methods to inte grat e the coupled-channel Schrodinger equation numerically. In this work we 
shall employ two packages MATSCS[42|| and FESSDE2.2[43 to perform the numerical calculation so that the results 
obtained by one program can be checked by another. The first package is a Matlab software, and the second is 
written in Fortran??. Both packages can fastly and accurately solve the eigenvalue problem for systems of coupled 
Schrodinger equations, and the results obtained by two codes are the same within error. 



The masses of the involved mesons are taken from PDG[4||: Md = 1869. 3MeV, Md* = 2006. ?MeV, Mdi = 
2422MeV, Mqo = 2308MeV, = 135MeV, = 54?.5MeV, = 600MeV, rup = ??5.5MeV and = 782.65MeV. 
The effective coupling constants in HM^PT have been studied from various phenomenological and theoretical ap- 
proaches, and the estimates for g, h are listed in Table |T1 It is obvious that there are still large uncertainties in their 
values. In the following, we shall first consider whether one pseudoscalar tt and 77 exchange can results in a bound 
state in the DiD and DqD* system, then the contribution of a exchange is included. 

The numerical results with only one pseudoscalar exchange are presented in Table |TT1 For several typical values of 
gh, we vary the cutoff A from a small value until we find a solution which lies below the DiD threshold. Here the 



TABLE I: Summary of theoretical estimates for the effective couphng g and h. 
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Reference 


g 


Remark 


[45j 


0.59 ±0.07 ±0.01 


combining the CLEO's results on D* decay width 


[46j 


0.46 ± 0.04 


through a constituent quark-meson model 


[47] 


0.53 


including one loop corrections without positive parity states 


[47] 


0.65 


including one loop corrections with positive parity states 


[27] 


0.44 ±0.16 


from QCD sum rule 


[48] 


0.39 ±0.16 


from QCD sum rule 


[49] 


0.32 ±0.02 




[28] 


0.75 


from non-relativistic quark model 



Reference h Remark 



[28] 


\h\ = 0.87 


from non-relativistic quark 


[50j 


0.9ll°o1 


in a constituent quark-meson model in soft pion limit 



mass M is measured with respeet to the DiD threshold Mp^ + ~ 4291.3 MeV, r^-n^s is the root of mean square 
radius, and R denotes the ratio between the DiD and DqD* components in the bound state solutions. By comparing 
the results with different A for the same value of the parameter gh, one notes that the magnitude of M increases with 
A, whereas the reverse is true for imis a-nd R. The bound state mass is sensitive to the parameter gh as well, larger gh 
is helpful to form a molecular state. From the numerical results in Table |TT1 we see that one can get a molecular state 
consistent with Y(4260), given appropriate value for gh and a reasonable cutoff A in the range 1-2 GeV. However, the 
existence of a bound state with I'^(J^) = l^(l^) require that the value of A should be at least larger than 4 GeV. The 
cutoff parameter A is a typical hadronic scale, which is generally expected to be in the range 1-2 GeV. If A is required 
to be much larger than 2 GeV in order to form a bound state, we tend to conclude that such a bound state should 
not exist. Therefore, it is not appropriate to assign Zj(4250) as a DiD or DqD* molecule, if only the non-diagonal 
interaction from tt and ry exchange is considered. 

Then we include the contribution coming from a exchange, which leads to only the diagonal interaction. The 
corresponding numerical results are shown in Table IIIII and Table IIVI The radial wavefunctions x('') = rR(r) for 
certain certain parameter values are shown in FiglS) The wavefunction corresponding to other solutions in Table 
IIIII and ITVl has similar shape with that in FiglJl We find that the a exchange interaction has significant effects, the 
variations of M, rrms and R with respect to A have the same pattern as those in the only pseudoscalar exchange case. 
Varying the parameters gh, gad'aj 9^9 a and h„ in the reasonable range results in large change of the predictions, 
which indicates that the results are sensitive to the effective coupling constants. We can see that large gh, negative 
gag'r, and 5^5^' are favorable to binding the molecular states. Comparing the results in Table [III IIIII and IIVI we find 
that a exchange is against the formation of bound state with I'^(J^'^) = 0~(1~~), nevertheless, it is beneficial to 
the formation of I^(J^) = 1~(1~) molecular state. As for Y(4260), the conclusion reached with only pseudoscalar 
exchange remains. Y(4260) could be accommodated as a molecule state for appropriate effective coupling constants 
and cutoff. A I'^(J^) = 1^(1~) bound state around 4250 MeV requires A should be at least 3 GeV, therefore we 
conclude that the interpretation of Z J (4250) as a DiD or DqD* molecule is disfavored. This conclusion is consistent 
with the general observations from chiral quark model. It is found that the isoscalar channel is easier to bind that 
the isovector channel for the same components [sij . 



B. The bottom analog of Y(4260) and Z^(4250) 

The bottom analog Ybb and Zj^j^ respectively denote the states obtained by replacing both the charm quark and 
antiquark with bottom quark and antiquark in Y(4260) and Zj(4250). The above calculation can be easily extended 
to study these states. The shape of both the diagonal and non-diagonal potential is similar to that of the charm 
system, except that the former is larger than the latter in magnitude. Furthermore, Since the kinetic energy is greatly 
reduced because of the heavier mass of B meson, a molecular state is more easily formed. We choose the same set 
of parameters as in the previous section. The numerical results with only pseudoscalar tt, 77 exchange are shown in 
Table IVl and the results with both pseudoscalar and a exchange are listed in Table IVTl and IVlIl As is expected, the 
magnitude M of the bottom analog is larger than that of the corresponding charmed state for the same parameters. 
The variation of M, rj-ms and R with A is the same as the charm system, large gh, negative g^g'^ and gada are beneficial 
to molecule formation as well. From the results in Table |Vl IVII and IVIII we note that both the bottom analog Ytb 
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r(fm) r(fm) 

(a) (b) 

FIG. 3: (color online) The radial wavefunction x('') = rR{r) for the molecular states with I*^(J^'^) = 0^ (1 ) and 
l''(J^) = l^(l^) respectively, (a) corresponds to the former state, and (b) is for the latter. We have taken gh — 0.85, 
gag'a = 0.58, go-flCT = 0.58 and \ha \ = h'„ = 0.35; A is chosen to be 1.4 GeV and 3.4 GeV respectively. 



and Zj^i^ may exist. 

Since Y(4260) has a large branch ratio into J / ip , the bottom analog Ybb should be searched for in the tt+tt^T 

channel. Specifically, the state Ybb can be searched for at B factories and future Super B factory via initial state 
radiation e^e~ 'yisR'^~^T^~'^ or by e"*"e~ — > tt+tt^T direct scanjs^. Furthermore, Ybb may be searched for at 
Tevatron via pp Ybb tt+tt^T, and LHC is more promising. Similarly, for the bottom analog Zt^, the most 
hopeful discovery channel would be "^^Xbi, where Xbi is in turn detected by its decay into 7T[4^. Because 

of the large mass of this state, it is difficult to produce such state via decay of certain particle ( i.e., Z J (4250) is 
produced in B decay ), consequently large hadron collides such as Tevatron and LHC are good place to search for 
this state. 



VI. CONCLUSION AND DISCUSSIONS 



In this work, we have performed a dynamical study of Y(4260) and Z J (4250) simultaneously to see whether they 
could be DiD or DqD* hadronic molecule. We have employed the HM^PT, which combines the heavy quark symmetry 
and the chiral symmetry. Since both the heavy meson and heavy anti-meson arc involved, the interaction related 
with heavy anti-meson has been included explicitly, and the total effective Lagrangian is invariant under the charge 
conjugation transformation. 

The off-diagonal interaction from pscudoscalar tt, rj exchange plays a dominant role, which is a straightforward 
support to the off-diagonal interaction mechanism proposed by Swanson and Close, a exchange leads to only diag- 
onal interaction, its contribution has been taken into account in this work. We find that a exchange is not favor- 
able to the formation of molecular state with I^(J^^) = 0~(1 ), whereas it is helpful to the binding of molecule 
with I^(J^) = 1~(1~). For appropriate value of the effective coupling constants and a reasonable cutoff A, Y(4260) 
could be accommodated as a DiD and DqD* molecule. However, the existence of a molecule around 4250 MeV 
with I'^(J^) = 1^(1^ ) requires that A should be at least 3 GcV, given the currently allowed values of the coupling 
constants. Consequently, the interpretation of Z^(4250) as a DiD or DqD* molecule is disfavored. Its structure 
should be studied further. Through calculating the masses of excited heavy tetraquarks with hidden charm in the 
diquark-antidiquark picture, the authors in Ref.fH^] suggested that Z J (4250) could be the charged partner of the 1^ 
IP state SS or as the 0~ IP state of the (SA ± SA)/-\/2 tetraquark. QCD sum rule analysis for Z J (4250) is performed 
in the Ref.Q- 

The effective potential from vector meson p, uj exchange has been presented analytically. Because of the accidental 
coincidence of mp and rriuj, the contribution from p and uj exchange almost cancels in the potential related with 
Z^(4250). For Y(4260), the situation is not the same. A number of effective coupling constants are involved. Because 
some of them have not been determined so far, we can not give a quantitative estimate about the contribution from 
vector meson exchange. Qualitatively, it should be smaller than the contribution coming from pscudoscalar and a 
exchange in magnitude. It is necessary and interesting to examine the effect of vector meson exchange on Y(4260) in 
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future. 

The bottom analog of Y(4260) and Z J (4250) denoted by Ybb and Z^^^ respectively may exist. Ybb can be searched 
for in e+e^ 7/sfi7r'''7r~T or by e+e^ tt+tt^T direct scan. The direct production of Ybb at Tevatron or LHC via 
pp Ybb ^ TT+TT^T is a hopeful approach as well. For Z^j^, the most promising discovery channel is Z^^ T^^Xbi- 
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APPENDIX A: THE POTENTIAL FROM LIGHT VECTOR MESONS p AND uj EXCHANGE 



The light vector mesons nonet can be introduced by using the hidden gauge symmetry approach, and the Lagrangian 
containing these particles is as follows [27l [3^. [33| 

+Mi(Tf )^''^^^.(p),,i?f ^) + h.c] (Al) 
where F^i,{p) = d^p^ - d^p^ + [p^, p^], and is defined as 

= ^^V, (A2) 
is a hermitian 3x3 matrix analogous to Eq.([5]) containing p, K* , lo and 0, 



V 



P^ 



V2 ^ V2 

K 




(A3) 



By imposing the KSRF relations, one obtains gv — 5.8. For the same reason, the interaction between the light vector 
resonances and heavy anti-mesons should be included via applying charge conjugation transformation 



+p,{H!^\-^F,ApUTp'')+h.c. 



(A4) 
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where we have used the property CV^C~^ = ■ Then the effective interactions reievant to the concerned tree level 
scattering diagrams are as follows 

+9u'DoV^laV^.abD7 + 9h'T5ov(^b'^ ^'^ ^ 9" DL)(a^K - d.V^)ab + /i.e. (A5) 
The coupling constants are as follows [ssj 



1 

ffD-DoV = .9d*DoV = -C5vv/2Md*Mdo 

1 

3 



3D-D„y = .9D-D„y = -"7^ (A6) 
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From the above effective interactions, following the general procedure presented in section III, we can calculate the 
effective potential from p and uj exchange. For Y(4260), the potential in coordinate space is 



vrr{r) = v^r{r) = o 



32M2,M2,^ 



3G(A,mp,r) + G{A,m^,r) 



1 2 



Mi(M2,-M2,J 



3F(A, X2,H5, mp, r) + Y{k, X2, ^e, rn^, r) 



+ 



ff^/xf(MD + MDj2 



72MdMdi 



2^MdMdi 

3Z(A, X2, /i5, TOp, r) + Z(A, X2, /Lt6, m^, r) 
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1 



^Z(A, X2, /is, TOp, r) H 5-.^ (A, X2, Ai6, "^c^, J') 



^22" W = [3iJ(A,mp,r) +i7(A,m„,r) 



32M2,.M2, 



3G(A, TOp, r) + G(A, r) 



1 2 
-^9v 



^Md-Mdo 



3F(A, X3, /i7, mp, r) + Y{A, X3, p,%,ma,, r) 



12Md.Md„ 

gie r 3 



3Z(A, ^3, fj.7,mp, r) + Z(A, X3, p,s, m^, r) 



12 



Z{K, X3, /i7, TOp, r) H 5-^(A, X3, /is, m^, r) 



(A7) 



The potential in coordinate space for Z J (4250) is 

V,T{r) = -i/3/325y[i?(A,TOp,r) -i7(A,c.,r) 



G(A, ?7ip, r) - G(A, ■m^,r) 



1 2 
-q9v 



Ci + 



/ii(M2,-M2,J 



2VMdMdi 



y(A,X2,/i5, Wp,r) - y(A, A:2,Ai6, w^^,?") 



72MdMdi 



Z(A, X2,/i5,TOp,r) - Z{A,X2,iJ.6,rnu;,r) 
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1 



Z(A, ATa, p.5,mp, r) 5-^(A, ATs, Me, 'ti^, 



^2^2" (^) = - [if (A, TOp, r) - (A, to^, r) 



32M2,.M2 



G(A, nip, r) - G(A, m^,r) 



L 2 
75y 



/i(M2). - M' 



2 1 

Do J 



2 r 



y(A, X3, /iy, "^p, - 5^(A, A:3, Ms, 'tiw, f ) 



^Md.Mdo 

Z(A, A:3, /i7, TOp, 7') - Z(A, A:3, ^s, m^-, r) 



12Md.Mdo 

T^(A, X3, M7, mp, r) T^(A, ^3, ^8, "i^, r) 

12 Lm; 



(A8) 



where the parameters p,i{i — 5, 6, 7, 8) are given by 



M^ = to2 - (Md, - Md)^ 

/i2 = to2 - (Md, - Md)2 

M? = to2-(Md„ -Md.)' 
M^=to2 -(Mdo-MdO' 



(A9) 
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The new functions G{A, m, r) and Y{A, X, /i, m, r) arc defined as follows 

G(A,m,r) ^ — l(A^e-^'' - m^e-"^'') + ' (A - -je'^'' 
47r r Stt r 

y(A,X,^,m,r) = i-i(e---e-^'-)- ^^i-^e-^^ (AlO) 

As is demonstrated in Eq. (jA7|) and Ea. (jA8[) . light vector mesons p and lo exchange leads to diagonal interaction, and 
the off-diagonal components of the effective potential are zero because of parity conservation. For Z J (4250), it is 
obvious that the potential coming from p exchange almost cancels that from oj exchange, because of the accidental 
coincidence of rup and m^, i.e., rup ~ 775.5 MeV and ~ 782.7 MeVQ. 

There are a number of parameters /3, /3i, (32^ P-i, C ^'^d Ci involved in the effective potential. The information 
about the effective coupling constants between the heavy meson and the light vector mesons is very scarce until now, 
especially those related with the P-wave heavy mesons. By vector meson dominance, /3 is estimated to be about 0.9 
[1^. Ref.[27t gives p, = — O.lGeV"^ and C, = 0.1. The remaining parameters have not been determined as far as we 
know, and we even don't know the ranges which they are in. So at present we can not give a quantitative estimate 
about the vector meson exchange contribution to the potential associated with Y(4260). Since the light vector meson 
mass TTip, m^^ is larger than m7r, and m„, we expect that the potential induced by vector meson exchange should 
be smaller than that due to pseudoscalar and scalar exchange in magnitude. In principle, we can determine these 
coupling constants following the methods of QCD sum rule, non-relativistic potential model and so on, by means of 
which certain coupling constants in HM^PT have been estimated. In future, if we could get a reliable estimate about 
these coupling constants from both phenomenological and theoretical approaches. The effective potential arising from 
p, LO exchange and its effect on the structure of Y(4260) could be analyzed in the same way as in section V. 
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1 cc 


with IG(JPC 


) = (r 
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Z+ with 


= 1 (r 


') 


gh 


A(GeV) M(MeV) 


rrms(fm) 


R 


A(GeV) M(MeV) 


rrms(fm) 


R 




3.3 


-4.04 


1.39 


2.63 


14.5 


-2.40 


1.76 


3.27 


0.23 


3.4 


-12.20 


0.82 


1.71 


14.6 


-7.53 


0.98 


2.02 




3.5 


-23.86 


0.63 


1.41 


14.7 


-14.71 


0.71 


1.60 




2.3 


-3.29 


1.56 


2.90 


9.6 


-3.53 


1.45 


2.76 


0.35 


2.4 


-11.32 


0.86 


1.74 


9.7 


-9.24 


0.89 


1.87 




2.5 


-24.95 


0.63 


1.40 


9.8 


-16.99 


0.68 


1.54 




1.6 


-1.79 


2.17 


3.88 


6.3 


-5.53 


1.17 


2.29 


A /I 

U.04 


1.7 


-10.33 


U.ao 




6.4 


-12.13 


U.oU 


1 71 
i. < i 




1.8 


-24.66 


0.67 


1.41 


6.5 


-20.82 


0.63 


1.45 




1.2 


-7.43 


1.15 


2.11 


4.0 


-3.45 


1.50 


2.81 


0.85 


1.3 


-22.69 


0.76 


1.46 


4.1 


-9.28 


0.93 


1.88 




1.4 


-46.18 


0.57 


1.24 


4.2 


-17.42 


0.71 


1.53 



TABLE II: The mass, the root of mean square radius(rms) and the ratio(R) between the DDi and D*Do components for the 
bound state solutions of the DDi and D*Do system with one pseudoscalar exchange, and the mass is measured with respect to 
the DiD threshold Md + Mdj ~ 4291. 3MeV. 
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1.88 


4 1 


^1 1 4 

-O J L4 


0.53 


0.90 


0.23 








4.6 


-7.02 


1.02 


1.80 


5.9 


-8.06 


0.87 


0.95 








0.35 


4 7 


1 4n 


0.72 


1.41 


fi n 

u.u 


1 8 ^9 


0.61 


0.73 




-0.58 


0.58 




4 S 


9fi "^4 

^U.04: 


0.58 


1.25 


fi 1 


"^1 "^4 

O -L .04 


0.50 


0.64 










8.4 


-3.81 


1.37 


2.32 


3.9 


-11.21 


0.74 


0.65 








0.50 


o.u 


q "^7 
y .o i 


0.88 


1.67 


4 n 


9^ 48 

.^0.40 


0.54 


0.50 










o.u 


1 fi 74 

±U. / 4: 


0.68 


1.42 


4.1 


4"^ 7fi 

40. 1 U 


0.45 


0.43 










4.4 


-6.47 


1.11 


2.32 


5.7 


-5.35 


1.15 


1.69 








0.35 


4.5 


-13.31 


0.80 


1.78 


5.8 


-12.98 


0.75 


1.17 




-0.58 


-0.58 




4.6 


-22.33 


0.64 


1.53 


5.9 


-23.10 


0.59 


0.96 










8.0 


-3.13 


1.60 


3.37 


3.8 


-8.99 


0.88 


1.12 








0.50 


8.1 


-7.05 


1.06 


2.36 


3.9 


-20.74 


0.62 


0.78 










8.2 


-12.35 


0.82 


1.92 


4.0 


-36.29 


0.50 


0.64 










2.8 


-2.15 


1.91 


3.38 


5.2 


-8.56 


0.90 


1.48 








0.35 


9 q 


q 4S 


0.93 


1.85 


o.o 


1 q 71 
-Ly . 1 ± 


0.62 


1.13 




0.58 


0.58 
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^U. ( u 


0.66 


1.46 


^ 4 

0.4: 


-34.24 


0.50 


0.98 










3.7 


-6.18 


1.12 


2.19 


3.6 


-3.77 


1.34 


1.87 








0.50 
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0.76 


1.64 
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0.72 


1.07 










q 


9^1 Q1 


0.59 


1.42 


R 
o.o 


4fi 

OU.4U 


0.53 


0.85 










2.8 


-9.66 


0.95 


2.15 


5.0 


-6.18 


1.12 


2.38 








0.35 


9 q 


1 q S9 


0.70 


1.70 


^1 1 

0. 1 


1 4 49 


0.76 


1.68 




0.58 


-0.58 




o.u 


-oo.^o 


0.56 


1.48 


9 


9^ ^R 


0.58 


1.38 










3.5 


-3.96 


1.45 


3.20 


3.5 


-4.65 


1.29 


2.55 








0.50 


o.u 


1 n Dfi 


0.93 


2.18 


fi 

o.u 


1 q7 


0.77 


1.54 










7 

O. 1 


1 S fifi 


0.70 


1.77 


7 

O. 1 


97 47 

-Zi I .4 ( 


0.57 


1.18 


0.35 








2.8 


-8.07 


0.99 


1.74 


5.0 


-7.24 


0.95 


1.25 








0.35 


9 q 


1 S 91 


0.70 


1.36 


^1 1 

0. 1 


1 fi q4 
- iu.y4 


0.66 


0.94 




-0.58 


0.58 




o.u 


■^1 84 


0.56 


1.21 


^ 9 


9q "il 
-zy .oi 


0.53 


0.82 










3.6 


-7.83 


0.99 


1.77 


3.5 


-5.98 


1.03 


1.19 








0.50 


3.7 


-16.23 


0.72 


1.42 


3.6 


-17.51 


0.65 


0.80 










3.8 


-27.23 


0.59 


1.26 


3.7 


-33.14 


0.51 


0.66 










2.7 


-8.08 


1.03 


2.05 


4.9 


-11.84 


0.81 


1.44 








0.35 


2.8 


-17.20 


0.74 


1.59 


5.0 


-21.48 


0.63 


1.17 




-0.58 


-0.58 




2.9 


-29.48 


0.59 


1.38 


5.1 


-33.64 


0.52 


1.03 










3.4 


-5.08 


1.27 


2.51 


3.4 


-5.68 


1.14 


1.80 








0.50 


3.5 


-11.34 


0.87 


1.84 


3.5 


-15.32 


0.72 


1.16 










3.6 


-19.83 


0.69 


1.55 


3.6 


-28.71 


0.56 


0.93 



TABLE III: The mass, the root of mean square radius(rms) and the ratio(R) between the DDi and D*Do components for the 
bound state solutions of the DDi and D*Do system with both one pseudoscalar exchange and a exchange, and the mass is 
measured with respect to the DiD threshold Md -I- Mdj ~ 4291. 3MeV. 







with IG(JPC) 


= (r 


) 


Z+ with I'^(JP 


) = 1 (1 




gh 




9a -da 


ha 


A(GeV) M(MeV) r 


rms(fm) 


R 


A(GeV) M(MeV) 


rrms(fni) 


R 










1.8 


-2.42 


1.84 


3.21 


4.2 


-9.99 


0.86 


1.59 








0.35 


1 q 


1 1 0"? 


0.91 


1.75 


4 % 

4.0 
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0.62 
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0.58 


0.58 




9 n 
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0.66 


1.39 
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0.50 
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2.1 
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3.2 


-12.62 


0.78 
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0.50 
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0.71 
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0.57 


1.08 










9 
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0.57 
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4 
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0.47 
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-6.29 
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2.42 


4.1 
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0.78 
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1 Q 
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0.78 
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94 8'^ 

-Z4.00 


0.61 
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0.58 
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0.60 
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4 
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^8 7q 
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0.50 


1.32 










2.0 


-5.65 


1.25 


2.51 


3.1 
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0.85 


1.87 








0.50 


9 1 


1 4 
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0.82 


1.77 


9 


94 no 


0.62 


1.43 










9 9 
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0.64 


1.49 


o.o 
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0.51 


1.21 
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-5.35 
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2.11 


4.1 
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0.74 
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0.78 
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0.60 
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0.79 
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0.50 
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1 4Q 


0.83 
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0.64 
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0.98 
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0.98 


1.84 
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-22.04 


0.70 


1.46 


4.0 


-16.85 


0.71 


1.42 




-0.58 


-0.58 




2.0 


-39.09 


0.56 


1.28 


4.1 


-28.01 


0.58 


1.22 










2.0 


-9.66 


0.98 


1.87 


3.0 


-10.72 


0.87 


1.60 








0.50 


2.1 


-20.31 


0.72 


1.48 


3.1 


-22.23 


0.64 


1.22 










2.2 


-34.80 


0.58 


1.31 


3.2 


-37.53 


0.52 


1.04 










1.3 


-11.86 


0.93 


1.72 


3.2 


-10.84 


0.86 


1.67 








0.35 
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98 TK 


0.67 


1.34 


o.o 


-22.29 


0.64 


1.34 




0.58 
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1 
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1 4 
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1.52 
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0.65 


1.32 
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0.52 


1.13 
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-14.66 


0.87 


1.71 


3.1 


-11.85 


0.86 


1.99 








0.35 


1 4- 


■^9 74 


0.64 


1.38 


^ 9 
o.z 


99 "^n 


0.66 


1.61 




0.58 


-0.58 




1 


^8 fi^ 


0.52 


1.24 


o.o 


q8 
-oo.yo 


0.54 


1.41 










1.3 


-7.58 


1.14 


2.13 


2.5 


-7.56 


1.07 


2.38 








0.50 


1.4 


-20.52 


0.76 


1.53 


9 (S 


1 7 q4 

-L 1 .y4 


0.73 


1.70 










1 


■^Q 78 
-oy. f o 


0.60 


1.31 


9 7 


^9 4q 
-oz.4y 


0.57 


1.40 


0.85 








1.3 


-14.13 


0.87 


1.56 


3.1 


-10.85 


0.85 


1.44 








0.35 


1 4 
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0.64 


1.26 


^ 9 
o.z 
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-z± .uy 


0.64 


1.17 




-0.58 


0.58 




1 


^8 1 fi 


0.52 


1.13 


o.o 


8n 


0.53 


1.05 










1.3 


-7.04 


1.15 


1.95 


2.5 


-6.58 


1.08 


1.74 








0.50 


1.4 


-19.98 


0.76 


1.38 


2.6 


-17.47 


0.71 


1.23 










1.5 


-39.26 


0.59 


1.20 


2.7 


-32.64 


0.56 


1.03 










1.2 


-4.75 


1.40 


2.51 


3.0 


-11.33 


0.87 


1.75 








0.35 


1.3 


-16.97 


0.82 


1.56 


3.1 


-21.18 


0.67 


1.42 




-0.58 


-0.58 




1.4 


-36.26 


0.62 


1.29 


3.2 


-33.99 


0.56 


1.25 










1.3 


-9.47 


1.03 


1.87 


2.5 


-13.90 


0.81 


1.62 








0.50 


1.4 


-23.58 


0.72 


1.41 


2.6 


-26.68 


0.62 


1.30 










1.5 


-44.09 


0.58 


1.23 


2.7 


-43.54 


0.51 


1.14 



TABLE IV: The continuing of Table Hill 
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Ybb 


with 


) = (1 


) 


Zjb 


with 


) = 1 (1" 


") 


gh 


A(GeV) M(MeV) 


rrms(fm) 


R 


A(GeV) M(MeV) 


rrms(fm) 


R 




1.8 


-5.60 


0.81 


3.08 


6.2 


-7.90 


0.63 


2.53 


0.23 


1.9 


-14.96 


0.53 


2.03 


6.3 


-13.45 


0.49 


2.04 




2.0 


-28.27 


0.42 


1.63 


6.4 


-20.22 


0.40 


1.76 




1.4 


-8.76 


0.70 


2.57 


4.2 


-5.17 


0.88 


3.08 


0.35 


1.5 


-21.93 


0.49 


1.79 


4.3 


-10.31 


0.57 


2.29 




1.6 


-40.65 


0.39 


1.48 


4.4 


-16.95 


0.45 


1.89 




1.1 


-11.29 


0.66 


2.35 


2.9 


-6.55 


0.73 


2.81 


A /I 

U.04 


1.2 


-29.05 


U.4 / 


i.DO 


3.0 


-12.79 


U.o4 


O 1 Q 

z.io 




1.3 


-54.93 


0.38 


1.38 


3.1 


-20.87 


0.44 


1.77 




0.8 


-3.49 


1.13 


4.14 


2.0 


-7.38 


0.72 


2.72 


0.85 


0.9 


-18.98 


0.58 


1.94 


2.1 


-14.97 


0.53 


2.03 




1.0 


-46.04 


0.43 


1.46 


2.2 


-25.03 


0.44 


1.69 



TABLE V: The mass, the root of mean square radius(rms) and the ratio(R) between the BBi and B*Bo components for the 
bound state solutions of the BBi and B*Bo system with one pseudoscalar exchange, and the mass is measured with respect to 
the BiB threshold Mb + Mbi ~ 11004MeV.. 
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0.67 
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3.4 
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1.9 


-14.61 


0.54 


1.85 


2.1 


-21.54 


0.42 


1.02 










2.0 


-25.57 


0.43 


1.56 


2.2 


-39.85 


0.35 
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1.5 


-8.80 


0.71 


2.68 


2.4 


-6.03 


0.78 


3.09 
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1.6 


-19.49 


0.51 


1.94 
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0.54 
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0.42 
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0.80 


3.08 
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0.76 


3.03 
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2.23 
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0.50 
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1.9 


-23.92 


0.47 


1.84 


2.1 


-31.57 


0.40 


1.27 



TABLE VI: The mass, the root of mean square radius(rms) and the ratio(R) between the BBi and B*Bo components for the 
bound state solutions of the BBi and B*Bo system with both one pseudoscalar exchange and a exchange, and the mass is 
measured with respect to the BiB threshold Mb -I- Mbi ~ 11004MeV. 





Ybb 


with IG(JPC 


) = (1 


) 


Zjb 


with 


) = 1 (1 


") 


gh 




9a -da 


ha 
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0.47 


1.60 










1 A 


-OtJ. lO 


0.43 


1.55 


9 n 


^fi ni 


0.37 


1.24 










1.1 


-6.49 


0.84 


3.19 


2.1 


-8.66 


0.68 


3.21 








0.35 


1 9 


1 Q 9^ 


0.55 


2.03 


9 9 


1 7 ^14 


0.50 


2.29 




0.58 


-0.58 




1 


■HS 49 


0.43 


1.62 


9 


9Q "ifi 
zy .iju 


0.41 


1.83 










1.2 


-10.20 


0.69 


2.60 


1.7 


-5.05 


0.88 


4.32 








0.50 


1 


9'^ 1 9 


0.51 


1.90 


1 8 


1 4 "^fi 


0.56 


2.48 










1 4- 


41 94 


0.41 


1.60 


1 Q 
1 .y 


98 "^8 


0.43 


1.79 


0.54 








1.1 


-5.38 


0.88 


2.98 


2.1 


-5.49 


0.76 


2.23 








0.35 


1 9 


1 7 Q4 


0.55 


1.81 


9 9 


1 4 7fi 


0.50 


1.49 




-0.58 


0.58 




1 


■^7 1 1 


0.43 


1.44 


9 % 


-27.44 


0.40 


1.20 










1.2 


-8.84 


0.71 


2.35 


1.8 


-12.06 


0.56 


1.62 








0.50 


1 


91 fi'H 


0.51 


1.69 


1 Q 
-L .y 


97 09 
.i 1 .uz 


0.41 


1.16 










1.4 




0.41 


1.42 


9 n 


47 (Ti 


0.34 


0.96 










1.1 


-7.79 


0.78 


2.81 


2.0 


6.63 


0.76 


3.06 








0.35 


1.2 


-21.69 


0.52 


1.85 


2.1 


-14.74 


0.54 


2.07 




-0.58 


-0.58 




1.3 


-42.15 


0.42 


1.50 


2.2 


-25.89 


0.43 


1.62 










1.2 


-12.19 


0.65 


2.30 


1.7 


-8.44 


0.69 


2.75 








0.50 


1.3 


-26.32 


0.49 


1.73 


1.8 


-20.06 


0.49 


1.76 










1.4 


-45.74 


0.40 


1.49 


1.9 


-36.41 


0.39 


1.35 










0.9 


-15.10 


0.63 


2.09 


1.8 


-14.35 


0.54 


1.93 








0.35 


1 n 


■^7 


0.46 


1.53 


1 q 
-L .y 


97 7^1 


0.42 


1.53 




0.58 


0.58 




1.1 




0.37 


1.31 


9 n 




0.35 


1.32 










0.9 


-12.39 


0.67 


2.24 


1.5 


-5.80 


0.80 


2.90 








0.50 


1 n 


■^1 sn 

-o± .ou 


0.48 


1.60 


1 fi 


1 7 fi4 


0.51 


1.79 










1.1 


-60.44 


0.39 


1.36 


1 7 


flQ 
Otj.uy 


0.40 


1.39 










0.9 


-16.45 


0.62 


2.10 


1.7 


-11.83 


0.62 


2.67 








0.35 


1 n 


4n 1 Q 


0.45 


1.56 


1 8 


9"^ 08 


0.47 


2.00 




0.58 


-0.58 




1 1 

J L 


74 Ql 


0.37 


1.35 


1 Q 
1 .y 


■^8 1 
-00. 10 


0.39 


1.66 










0.9 


-13.70 


0.65 


2.25 


1.5 


-11.33 


0.64 


2.73 








0.50 


1 n 


■^4 DQ 


0.48 


1.64 


1 fi 


-24.92 


0.47 


1.91 










1 1 
J L 


-uo.oo 


0.39 


1.40 


1 7 


4"^ 8^ 

-4:0.00 


0.38 


1.53 


0.85 








0.9 


-15.93 


0.62 


2.01 


1.7 


-9.30 


0.64 


2.04 








0.35 


1 n 


-oy.ou 


0.45 


1.47 


1 8 


9n 77 


0.47 


1.49 




-0.58 


0.58 




1 1 

J L 


74 9^ 


0.37 


1.27 


1 Q 
1 .y 


■^fi 9n 


0.39 


1.25 










0.9 


-13.17 


0.66 


2.15 


1.5 


-9.18 


0.65 


2.11 








0.50 


1.0 


-33.43 


0.48 


1.54 


1.6 


-23.12 


0.46 


1.44 










1.1 


-63.13 


0.39 


1.31 


1.7 


-42.63 


0.38 


1.16 










0.9 


-17.27 


0.61 


2.02 


1.6 


-7.05 


0.77 


3.02 








0.35 


1.0 


-41.88 


0.45 


1.51 


1.7 


-16.42 


0.54 


2.04 




-0.58 


-0.58 




1.1 


-77.70 


0.37 


1.30 


1.8 


-29.47 


0.43 


1.61 










0.9 


-14.47 


0.64 


2.16 


1.5 


-14.85 


0.57 


2.14 








0.50 


1.0 


-35.71 


0.47 


1.57 


1.6 


-30.27 


0.44 


1.57 










1.1 


-66.52 


0.38 


1.35 


1.7 


-51.10 


0.36 


1.30 



TABLE VII: The continuing of Table IVll 



